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ABSTRACT 

Some features of the physics of radiation-dominated shock waves are discussed with emphasis 
on the peculiarities which are important for correct numerical modeling of shock breakouts 
in supernova. With account of those peculiarities, a number of models for different supernova 
types is constructed based on multigroup radiation transfer coupled to hydrodynamics. We 
describe the implementation of a new algorithm RADA, designed for modeling photon trans- 
fer at extremely-relativistic motions of matter, into our older code STELLA. The results of 
numerical simulations of light curves, and continuum spectra are presented. The influence of 
effects of photon scattering on electrons, of thermalization depth and of special relativity in 
transfer equation is considered. Some cases are demonstrated, when the appearance of hard 
X-ray emission is possible at the shock breakout. The necessary refinements in numerical 
algorithms for radiative transfer and hydrodynamics are pointed out. Prospects for using the 
results of numerical simulation to analyze and interpret available and future data from space 
observatories are discussed. 
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1 INTRODUCTION 

The phenomenon of a supernova in most cases should start with a 
bright flash, caused by a shock wave emerging on the surface of the 
star after the phase of collapse or thermonuclear explosion in inte- 
riors. The detection of such outbursts associated with the supernova 
shock breakout can be used to obtain information about the explo- 
sion properties and presupemova parameters, which is necessary to 
understand the physical processes that underlie this phenomenon. 

For an accurate treatment of the shock wave propagation near 
the surface of a presupemova it is necessary to perform numerical 
calculations which, in addition to hydrodynamics, should account 
coiTectly for radiative transfer in moving media. In some cases, e.g. 
in compact type Ib/c presupernovae, shock w aves can reach rela- 
tivistic velocities jBlinnikov et alj2002ll2003h . Then one has to in- 
clude into consideration a number of relativistic effects. This paper 
presents the results of numerical simulation of shock waves in sev- 
eral models for type lb and type II supernovae, taking into account 
these features. 

Studying the supernova shock breakout becomes particu- 
larly topical in connection with the recent detection of this phe- 
nomenon by the SWIFT spacecraft ( Soderberg __et al 20081). In ad- 
dition, a flash from the shock wave on t he surface of red super- 
giant in type II sup ernova SNLS-04D2dc jSchawinski et alj|200i ; 
iGezari etaL 20081) and light echo from the shock wave of Gas A 
(IPwek & Arendtll2008h are detected. Simulations jTominaga et al.l 
[2009.. 20 lH) with our code STELLA show good agreement with ob- 
servations in case of SNLS-04D2dc. The enormous luminosity of 



SN 2006gy is most successfully explained by models where the ra- 
diative shock wave provides almost all radiation for many months 
jWooslev. Blinnikov & Hegeill2007h . There is a possibility of ob- 
taining new data if the exp eriments similar to LOBSTER space 
observatory dCalzavara & M atzner 20q3), the satellite EXIST (En- 
ergetic X-ray Im aging Survey Telescope, see iGrindlav et al II2OO3I : 
iBand et alj[20o3) or any X-ray station of a similar type would be 
launched in the future. E.g., the experiment MAXI (Monitor of All- 
sky X-ray Image) on board the module Kibo at ISS jMatsuoka et al.l 
1199"^) is already started. 

In this paper we compare the results of the predictions for 
shock breakout for three codes: the equilibrium diffusion gray 
radiation-hydro code SNV, nonequilibrium multigroup radiation- 
hydro code STELLA and ultrarelativistic nonequilibrium multi- 
group transfer code RADA. We check the sensitivity of the pre- 
dictions of the shock emission to the parameters of the numeri- 
cal scheme, such as the boundaries of the frequency interval (in- 
cluding the X-ray range), and to the approximations in the opac- 
ity description. It is found that high-temperature peak behind the 
shock front and production of hard "tail" in spectrum are sup- 
pressed by extremely low true absorption, with a cross-section at 
the level 10~^ of Thomson scattering in a presupemova SN lb. 
This level of absorption can be provided by double Com pton effect 
jMandl & Skvrmelll952l : rWeaveJl976l:lLightmanlll98lh , so it must 
be taken into account in realistic models of radiation-dominated 
shock waves in a supemova. Some additional refinements neces- 
sary in the methods of constmcting shock-breakout models are 
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pointed out. Among them we consider how the relativistic and geo- 
metric effects in radiative transfer in a comoving frame of reference 
influence the predictions of supernova light curves and spectra at 
the epoch of shock-breakout. 



2 RADIATION DOMINATED SHOCK WAVES 

Let us consider some of the peculiarities of a shock wave when 
one cannot neglect radiation, and, moreover, when the shock front 
structure is determined by energy and momentum of photons. This 
situation is typical for supernova explosions, therefore, its descrip- 
tion is important for building adequate models. 

There are many theoretical articles and several books where 
properties of radiation-dominated shock are considered. In most 
theoretical works it is supposed that the medium is optically thick, 
i.e. photons bom in the shock waves are absorbed by the cold mat- 
ter upstream. This pictures holds in a stellar explosion before the 
actual outburst of supernova light while the shock is still buried 
under the photosphere. 

A brief review of the contri butions that are important fo r the 
history of the problem is given bv lBUnnikov & TolstovUioiTI) . Be- 
low we mention only the papers most important for further discus- 
sion. 

A very detaile d study of shoc k structure with radiation has 
been performed by IWeaveil (Il976l) : he has shown that produc- 
tion of observable hard (gamma-ray) radiation at supernova shock 
breakou t is very problematic. Nevertheless, Blandford & P avnd 
( Il981bl l3) have predicted the possibility of a nonthermal "tail" 
in spectra bom in shock waves in hot media. However, these 
authors have neglected a number of eff ects associated with 
Compton scattering taken into account by IWeavej (Il976h . Fur- 
ther refinement of the theory wi th regard to the thermal Comp - 
ton effect has been mad e by Lvubarskii & SvunvaevI ( Il982h . 
iFukue. Kato & Matsumot"3 tl985r ). iBeckej h988ll aiid other re- 
searchers. 

Discontinuities inevitably develop in gas parameters in a 
shock wave which is not too weak and not too strong (while ra- 
diation energy and momentum can be neglected), if one takes into 
account only thermal conduction and n eglects viscosity (see sec- 
tion VII. 3 in IZeldovich & Raizej[T963) . In real gases the discon- 
tinuities are smoothed by viscosity which transforms the ordered 
motion (kinetic energy of the flow) into the heat (chaotic motion 
of particles) on the distance of the orde r of one mean free path of 
particles. This is called a viscous jump dZeldovich & Raizeilll966l : 
[Shiilll992h . At very high amplitudes of the shock, when the energy 
density and pressure of radiation become large in comparison with 
those of matter, the situation is different and the viscous jump dis- 
appears. Gas goes from the initial state to the final one under the 
action of radiative heat conduction, even if the viscosity of mat- 
ter is not taken into account (Belokon 1959j). Such a situation is 
achieved when radiation pressure Pr and the gas pressure Pg he.- 



hind the shock front are rela t ed as P r / Pg — 4.4. 

Ilmshennik & Morozovl ( 1 19641) analysed various options for 
transfer equation (simple diffusion equation and the Eddington ap- 
proximation in the moving medium). Plasma upstream the shock 
front was taken to be cold (To = 0, Pg = 0). Their solution 
is more rigorous than presented by (lBelokonlll959b and it yields 



Pr/Po 



.5 for the condition of disappearance of viscous jump. 



The compression in the shoc k should then exceed a crit ical value 
Pi/po = 6.68. However, see lWeaver & Chaplind ( [l974l) . where it 
is again obtained Pr/Pg — 4.4. It should be noted that in determin- 



ing the stmcture of the shock front, the flux of photon momentum, 
i.e. the radiation pressure, begins to play a role comparable to the 
flux of energy. Such conditions are often achieved in astrophysical 
shock waves, in particular, in superaovae. 

The viscous density jump disappears during propagation of 
the shock wave inside of a presupemova if radiation pressure 
is appreciably higher than the plasma pressure and when there 
is a fairly intense processes of absorption and creation of pho- 
tons, not just their pure m onochromatic scattering. Note that 
Ilmshennik & Morozovl ( 1 1964b assumed extinction equal to the tme 
absorption, scatterin g was not taken into account. 

Calculations of lWeaveJ ( Il976l) have been performed in the 
diffusion approximation with a simplified description of Comp- 
ton scattering and various processes of absorption and creation 
of photons, and they have resulted in a critical ratio Pr/Pg — 
4.45. He assumed that the radiatio n has an equilibrium blackbo dy 
spectmm. Th i s was criticize d by iBlandford & Pavn3 l ll981bl lal): 
iRifferll ( ll988l) : lBeckeil ( Il988l) . In our work, we make no assump- 
tions about the equilibrium of radiation, and take into account in 
the calculations of the evo lution of photons the same terms as 
IBlandford & Pavnd (Il981bllal) Nevertheless, our results are closer 
to the results of lweaved ( ll976h . than to the results of his critics: we 
show that the achievement of high temperatures in the radiation- 
dominated shock waves is unlikely with a realistic allowance for 
the photon production processes. 

An important parameter is the optical depth of the 
heating zone for strong shocks. In the approximation of 
Ilmshennik & Morozovl ( 1 19641) , 



12c 



(1) 



where Ds is the shock front velocity and c is the speed of light. 
The approximate solution bv lKlimishirJjl968l) uses the concept of 
a critic al temperature Tfc , introduced by lRaizeil(ll957h . The solu- 
tion bv lKlimishirJ h968') is closer to the conditions of supemovae 
than the solution by Ilmshennik & Morozov C1964]) , because in 
supemovae is low at a low density. It gives 



14c 



(2) 



i.e., a value of the same order of magnitude. 

Thus, the shock wave propagates without energy losses 
as long as th e lower limit on the opt i cal depth holds 
(lOhvamal Il963l: I Imshennik & Morozovl Il964l: Ichevalieil Il98lt 
Ilmshennik & NadezhidI 19891) : 



- ^ > _L 



(3) 



where SR is the distance from the shock front to the photosphere, 
/ i s the photon m ean free path. Estimate ((S} was actually obtained 
bv lSachslll946l) as the thickness of the radiation-dominated shock 
front instead of the more accurate estimates |[T) and (|2). To obtain 
an estimate ([3}, the shock front travel time tsu = SR/Ds to the 
distance 5R should be set equal to td = iSR)'^ /cl — the photon 
diffusion time to the same distance. 

As the shock approaches the photosphere, relation (|3) breaks 
down and an outburst occurs at the stellar surface. The shock prop- 
agation in this regime can no longer be considered adiabatic, which 
makes it difficult to constmct analytical solutions and necessitates 
numerical calculations of a process in which radiative transport 
plays a very important role. 
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3 SELF-SIMILAR AND NUMERICAL SOLUTIONS OF 
SUPERNOVA SHOCK BREAKOUT 

3.1 Self-similar solutions 

The proble m of shock breakout in decreasing density was first for- 
mulated by iGandelman & Frank-Kamenetsk'vl jl95d) . They found 
a self-similar solution for a shock wave approaching stellar surface 
down the density decreasing as a power of distance to the surface. 
The properties of the solution for various values of the power and 
adiabatic index w ere calculated by iGrasberg (1981). The solution 
was confirmed bv lSakurail l ll960l) who showed in addition that sub- 
sequent expansion in a rarefaction wave can be described by an- 
other self-similar solution. The bulk of the pressure created by the 
shock before emergence accelerates matter in the rar efaction wave 
by a f actor of ~ (1.5 — 2) (see details in iLitvinova & Na dezhirt 
Ii99^). 

As the shock approaches the surface where the pre- 
supemova matter density p falls sharply, it accelerates fol- 
lowing a self-similar solution. A good overview of analyti- 
cal solutions and est i mates for this problem was given by 
llmshennik & Nadezhiij 119891) : newe r analytical results were pre - 
sented by 'Matzner & McKee' fl999^. 'Johnson & McKee| _tl97ll ): 
IColgate. M cKee & Blevins (1972); McKee & Colgate i 1973) stud- 
ied the shock acceleration to relativistic velocities in compact pre- 
supemovae. Previously, it was assumed that this could generate 
an outburst of X-ray or even g amma-ray radiation (IColgatelll974l ; 
ISisnovatvi-Kogan et alj 11973 ). The difficulty oi^the gamma-ray 
photon production in this scenario was shown bv iWeaver, U976I) . 



3.2 Numerical algorithm stella 

Very few detailed calculations of shock breakout have been pub- 
lished. Previously, they were made by invoking a number of rough 
approximations, such as the use of constant Eddington factors, 
the single-group ap proximation, and the neglect of the expan- 
sion effect in op acity ( Klein & Chevalier! 1978l : lEnsman & Burrows! 



I 991 I Kelly & Korevaai 199^). Here, we use the code STELLA 



I Blinnikov et alJll998Ll2006l) . which is designed to solve the prob- 



lem of the radiative transfer of nonequilibrium radiation with al- 
lowance made for hydrodynamics and, subsequently, to model the 
light curves of supemovae. STELLA is the hydrodynamics code that 
incorporates multigroup radiative transfer We have a veiy fine grid 
in outer layers and there are always several mesh points within the 
shock front. The absence of viscous jump in radiation-dominated 
shocks is vital for computing their hydrodynamics with accuracy 
of the computation of radiative transfer Thus, the use of artificial 
viscosity is not necessaiy for calculation of the radiation-dominated 
shocks in outer layers. 

The time-dependent equations are solved implicitly for the an- 
gular moments of intensity averaged over fixed frequency bands. 
The number of frequency groups available for current worksta- 
tion computing power, typically 100-300 in the range from lAto 
5 • 10* A, is adequate to represent, with reasonable accuracy, the 
nonequilibrium continuum radiation. STELLA includes in full opac- 
ity photoionization, free-free absorption, lines and electron scat- 
tering. The equation of state treats the ionization by equilibrium 
Sahas approximation. Thus, we use the method of complete multi- 
group radiation hydrodynamics in which the defects of older ap- 
proaches were corrected. This method is well applicable for the 
models of type II sup e mova e, such as SN 1987A and SN 1993J 
dBlinnikov etal]|l998l , I2OO0I) . Our STELLA method provides the 



most reliable predictions for an outburst to be made as long as 
the matter velocity u is less than ~ (20 — 30) % of the speed of 
light c. The method consistently takes into account all the terms of 
the order of u/c in hydrodynamic and radiation transfer equations 
whereas the terms of the order of {u/cf' are neglected. 



3.3 Numerical algorithm rada 

If the velocity of matter behind the shock wave reaches a significant 
fraction of the speed of light, j3 = u/c > (0.2 — 0.3), then it is 
more relevant for STELLA to use an algorithm RADA, which is able 
to solve the radiative transfer equation in comoving frame up to the 
values of the Lorentz factor, 7 = (1 — ~ 1000. The use 

of RADA is more relevant for relativistic flows in the most energetic 
supemovae and ORB afterglows, but in this paper we would like to 
use STELLA for testing RADA algorithm and find out the corrections 
to STELLA algorithm at the limit of app licability. Th e transport 



ty. Itie 
^[l980h : 



equation solved by RADA looks like this l lMihala! 

J(l + m'-^ + 7(. + P)'-^ + 7(1 - X 
(1 + l\.^a MR\^P 



+7W + /3)fj^ a. 

+ ^(l + /5M)f +7V(M + /3)f 
c Ot or 



(4) 



Here, r) - emission coefficient, x - absorption coefficient, /i - co- 
sine of the angle between the photon momentum and the radial di- 
rection. All values refer to the comoving frame of reference. We 
present here the equation in order to give the reader the scale of its 
complexity and non-triviality of its numerical solution. 

To simultaneously solve the relativistic radiative transfer equa- 
tion in a comoving frame of reference with the hydrodynamic equa- 
tions, we must determine the coupling, that is the quantities ap- 
pearing both in the transfer and in hydrodynamics. For example, 
the radiation moments Jv,Hv,Kv can act as these quantities. They 
should be found in the calculations at hydrodynamical grid points 
Ti, Vm, i.e., for each energy group Um and each Lagrangean zone 
ri. Thus, the problem is reduced to finding the radiation moments 
on the radius-frequency grid ri, Vm at each instant of time in the 
outer stellar layers where the optical depth is small. To calculate the 
moments, we must know the radiation intensity Ii^k,m as a function 
of fl. 

Here, we use a characteristic method of solving the radia- 
tive t ransfer equation developed previously (Tolstov & Blinniko^ 
I2OO3I) . This allows the radiation intensity to be found at a given 
grid point for some set of cosines /x;. The set of ^; is chosen by 
optimizing the calculation of the first moment Ji, with a specified 
accuracy by the number of points fii at a given grid point Vi.Vm- 

The initial intensities at the first step Ik are chosen under the 
assumption of blackbody radiation and at each successive step start- 
ing from tK+i are the result of a linear interpolation of the inten- 
sities calculated at the previous step. In outer stellar layers the op- 
tical depth is small, but the radiation becomes increasingly close to 
blackbody spectrum as one goes deeper into the star Therefore, the 
radiative transfer equation is numerically solved only in the outer 
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layers (as a rule, several tens of radial mesh points) and the radia- 
tion in deeper layers is assumed to be blackbody. 

For each grid point Vi, u,n, we solve the relativistic radiative 
transfer equation by the method of short characteristics. The char- 
acteristic emanates from point r^, jj., Prn, tk and ends at time tk+i, 
making some number of steps s„, n — 1...N. In our calculations, 
we assume a linear change in emission and absorption coefficients 
7] and X in the interval of the characteristic ds due to different times 
at the boundaries of this interval: 



dl{s) 
ds 



= ri + ris + xl{s) + xsl{s) , 



(5) 



where r],x,'ri,X constants. The analytical solution of this equa- 
tion is: 



/(.) = 



lo — — 

X 



+ 1 
X 



n 

X 



(6) 



where t — 0.5xs + X^- 

Apart from the radiative transfer equation in a comoving frame 
of reference, RADA takes into account the differences in delay of 
the radiation from the supernova explosion. This stems from the 
fact that the radiation from the edge of the star visible to an observer 
comes later than that from the central regions. 

Solving radiative transfer equation (|4) is important for al- 
lowance for time delay effect in the most energetic supemovae due 
to relativistic corrections. Allowance for this can not be performed 
by STELLA algorithm which solves radiative transfer by moment 
equations. 

The radiation flux at some time tabs at the point of ob- 
servation is an integral over all visible points of a stellar surface 
radiating with some intensity I determined by the solution of the 
transfer equation: 



F (f ^ - 2^ 



fiR^ {t) I[t, u, fi, R(t)] d^ . 



(7) 



Here, f is the radiation frequency, D is the distance of the star from 
a remote observer, /i is the cosine of the angle between the normal 
to the radiating surface and the line along observer's direction that 
intersects the stellar surface at local time t and at radius R(t). 
The observer's time tabs is connected with t by the relation 



D R 

tobe ~ t -\ fl . 

C C 



(8) 



Equations J7l8t are written in a reference frame at rest. The quanti- 
ties /, u, fl therein are connected by standard transformations with 
their values provided by RADA in comoving frame. 

The light curve at the epoch of supernova shock breakout 
could be affected significantly in case of large and highly accel- 
erated envelopes (see details in Appendix B). We will demonstrate 
examples of this in the next sections. 



4 SHOCK BREAKOUT IN SUPERNOVA MODELS 

We have already presented some of the results on the first outburst 
of hard radiation generated by the shock breaking out from the pre- 
supemova for the models of SN 1993J dSlinnikov et al.lll998h and 
for SN lb ( iBlinnikov et"aLll200l l2003h . Some of the results were 
only reported at conferences by S.I. Blinnikov in 

199iB 

Below, we 

publish these results and provide additional information about the 



^ |http://online.kitp. ucsb.edu/online/supemova/snovaetrans .html I 
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Figure 1. The model SNII 11.06 Mq. 



outbursts at shock breakout in SN II-P (in particular, the SN 1987 A 
type) and SN lb and consider the dependence of outburst properties 
and shock breakout hydrodynamics on presupemova parameters. 



4.1 Shock wave breakout calculated by SNV code 

Some properties of the shock wave breakout as calculated by 
D.K. Nadyozhin in 1993 are gathered in Table [T] The calcula- 
tions were done in the approximation of the radiation heat con- 
ductivity (that is equilibrium diffusion of radiation) with the aid 
of the hydrodynamic code SNV that was us ed in previous work of 
[Orassberg et al. ( 1971) and llmshennik & Nadezhin ( 1989b. All the 
presupemova models are due to lWeaver & Wooslevl(ll993b . 

The data for six models from Table[T]were additionally worked 
up to take into account the time-delay spread of the light curves. 
The algorithm of the corresponding filtering procedure is described 
in Appendix A. The results are compiled in Table |2] and shown in 
Figures [T}[6] below. In all the figures the dashed lines correspond 
to the time-spread light curves while the solid ones show the case 
when the time-delay effect is neglected. In order to make the light 
curves discernible both near the peak and in their tail part, the hor- 
izontal time axes are shown in terms of log(t — to) where to is the 
time the SW takes to reach stellar surface. 

Figure|7] shows the spectral luminosities for the SN lb model 
at 5 different points of time (as usual the dashed curves corre- 
spond to the time-spread case). The solid curves are nothing else 
but the Planck distributions at given effective temperature T{t). 
The time t = 1.04s corresponds to the luminosity peak with 
T{t) ^ Tp = 5.1-10'^ K (TableQ and the spectrum attains its 
maximum al hiy = 2.82 kT — 1.24 keV. The dashed curves rep- 
resent the result of the superposition of the Planckian spectra with 
the temperature T{t) varying in the time "window" [t — R/c, t]. 
Consequently, the time-spread spectral luminosity keeps a good ad- 
mixture of the keV photons during At ~ R/c that is much longer 
than the intrinsic luminosity does (compare the dashed and solid 
curves att = 2.51s). With time, the difference between the dashed 
and solid curves disappears: att = 11.2 s they become virtually 
indistinguishable. Note that the total (time-integrated) spectral lu- 
minosity is exactly the same for both the cases. 

In Figures [814101 the results are presented calculated by 
STELLA and RADA algorithms where radiation transfer is calcu- 
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Table 1. The intrinsic properties of the light curves during SW breakout. 



Model 


M 


R 


-^cxp 




Aip3 


hp 




(M0) 


(Rq) 


(10^1 erg) 


(lO^K) 


(sec) 




SN II 


11.06 


335 


1 3 


2 94 


1000 


2 92 


SNII 


15.08 


496 


1.3 


2.45 


1100 


3.08 


SNII 


15.08 


496 


1.8 


2.70 


720 


5.01 


SNII 


25.14 


881 


1.3 


1.53 


7200 


1.48 


SNII 


35.19 


1160 


1.3 


1.31 


6900 


1.37 


SN 1993 J 
(jl3a7) 


3.81 


629 


1.3 


2.70 


1900 


7.24 


SN lb 
trancated 
at zone 327 


3.51211 


0.763 


1.3 


51.0 


0.12 


1.45 


SN lb 
trancated 
at zone 346 


3.51250 


1.23 


1.3 


35.0 


0.60 


0.807 


SN 1987A 


22.12 


66.6 


1.3 


6.01 


90 


2.07 


SN 1987A 


18.10 


37.8 


1.3 


7.47 


33 


1.61 



M and R — presupemova mass and radius, respectively 



i?cxp — the kinetic energy at infinity 
Tp — the peak temperature 

Afp3 — the width of the light curve peak at T = i Tp (3 stellar magnitudes below 
the peak luminosity) 

Lp — the peak luminosity (in lO'^^erg s~^) 



Table 2. The comparison of the intrinsic and time-spread light curves. 



Model 


M 


R/c 




Lps 


Atpi 


Atpis 




(Mq) 


(sec) 


(lo-is 


erg s^^) 


(sec) 


(sec) 


SNH 


11.06 


779 


2.92 


1.09 


200 


672 


SNII 


15.08 


1150 


3.08 


1.27 


343 


991 


SN 11 


35.19 


2690 


1.37 


1.07 


2390 


3230 


SN 1993J 


3.81 


1460 


7.24 


3.20 


462 


1380 


SN lb 


3.51211 


1.77 


1.45 


0.056 


0.028 


1.42 


SN 1987A 


18.10 


87.8 


1.61 


0.358 


12.2 


59.4 



Here Lp and Lps are the peak intiinsic and time-spread luminosities, Atpi 

and Atpis are the widths of the light curves at 1 stellar magnitude be- 
low Lp and Lpg, respectively. For all the supernova models in Table |2] 
Sexp = 1.3-1051 erg. 

lated numerically. As in all the figures above the dashed lines cor- 
respond to the time-spread light curves while the solid ones show 
the case when the time-delay effect is neglected. The comparison 
with Figures |5}l7] shows that exact calculation of radiation transfer 
retains the qualitative form of the the light curves and spectra and 
adds a number of details. 



4.2 Shock breakout in the models of SN 1987A 

For our calculations of the SN 1987A outburst, w e used the presu- 
pemova model constructed by the Tokyo group dShigevama et al.l 
Il987l : IShigevama & Noniotolll990l). For details of this model and 
the va riants of our runs, see iBlinnikovl (Il999l) : | Blinnikov et al.l 
The variants of this series are designated as 14E1, 14E1.3 



etc., where the number 14 corresponds to the ejecta mass in solar 
masses A/0(to be more precise, 14.7 Mq). The numbers after E 
denote the explosion energy in units of 10^^ erg. 

Figures [TT] - [12] show the velocity, density and temperature 
profiles for one of these variants, 14E1.3. 

In the outermost layers, condition ([3) is violated, the losses 
through radiation become significant, and the shock acceleration 
predicted by the self-similar solution ceases. This behavior is ex- 
cellently seen from Fig. [13] where the velocity profiles are shown 
as a function of Lagr angean mass and optical depth r in the model 
I4E1 for SN 1987A ( lBUnnikov|[T999l) . The "breakout" of photons 
from behind the shock front occurs precisely at r ~ 10. These 
photons slightly accelerate the overlying layers of matter, but the 
cumulation of energy at a low mass is already inefficient due to 
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Figure 2. The model SNII 15.08 Mq. 
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Figure 5. The model SNlb 3.51211 Mq. 
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Figure 3. The model SNII 35.19 Mq. 
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Figure 4. The model SN 1993J 3.81 Mq. 
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Iog(t-5.l45x10^s) 
Figure 6. The model SN I987A IS.IOMq. 

great losses through radiation. The asymptotic ejecta ve locity dis- 
tribution in mass for this variant is given in the table from lBlinnikovl 
(fT999,). 

The velocity distribution of the ejected mass, the mass spec- 
trum, is of great importance in astrophysical applications. The mass 
spectrum M(u) with a velocity greater than u can be predicted an- 
alytically at the self-similar stage. Analysis shows that the mass 
spectrum for the outer part of the envelope is defined by a simple 
relation. It obeys the Nadyozhin-Frank-Kamenetskii law 



M(u) = const • u ^ 



(9) 



where M (it) is the mass with a velocity greater than u, the numer- 
ical va lue of the exponent is written out for t he case of interest to 
us fsee lNadezhin & Frank-Kamenetskiil965h . n is the exponent in 
the law p (X {R — r) " l|Grasberg"l981^ and A is the e igenvalue of 
the self-simil ar problem (Gandelman & Frank-Kamen etskv 195^; 
ISakurailll960h . 

The mass of the outer part of the envelope described by law 
(|9j is determined by the applicability conditions for self- similar so- 
lutions. For example, our calculation gives it = 5.5 x 10^ km/s at 
M = O.IMq. The prediction of u = 30 x lO^km/s must be at 



© 2012 RAS, MNRAS QOQ.miSol 



Coupling of matter and radiation at supernova shock breakout 7 




14ElshbaXrada 




lg(t-6.55xl03s) 
Figure 8. The model SN 1987A 16.3 Mq 



Figure [TT] shows similar velocity profiles for the variant 
14E1.3, where the explosion energy is higher than that in 14E1 
by 30%. The picture is very similar to Fig. [13] but now the break- 
out time decreased and the asymptotic ejecta velocity slightly in- 
creased. It is interesting to see how this picture appears in Eule- 
rian coordinates (Fig. I lit. Figure [T4lshows how the density peak is 
formed in the outer layers due to inefficient gas acceleration. Our 
numerical experiments show that the development of one, two, or 
more peaks is possible: this is analogous of the loss of stability of 
sequential harmonics. In this case, we see them at the nonlinear 
stage. How many peaks there are and how strong they are depends 
on the shock strength, on the density profile, on the absorption and 




hi^, keV 



10. Spectral luminosities for the model SNIb/c (see Figj2]for the 
light curve) at time stamps related to the maximum of the flux. 

emission coefficients, and (in our calculations) on the grid fragmen- 
tation and artificial dissipation parameters. This question could be 
investigated, but such a problem would be purely academic in the 
one-dimensional approximation — because even the formation of 
the first peak must lead to the development of multidimensional 
instabilities in it (in a thin layer) and to layer fragmentation. An- 
other analogy is a steepening of a great wave coming on shore; it 
can be seen how increasingly high harmonics are formed from one 
sine wave followed by breaking — everything flies away in small 
splashes (multidimensionality). 

Let us now consider the influence of opacity on the parame- 
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Figure 11. Matter velocity for the valiant 14E1.3 at the epoch of shock 
breakout versus Lagrangean ma ss Mr (top ) and E ulerian radius r (bottom) 
in the model for SN 1987A from lBlinnikovl j 19991) . The proper time is given 
near the curves. 



Figure 12. Matter temperature for the versions 14E1.3 (top) and 14E1X2 
(bottom) at shock breakout versus Lagrangean mass Mr measured from the 
surface. The time in seconds is given near the curves. 



ters of the shock breaking out from a presupemova. We see from 
FigJ12lthat the matter te mperature behind the s hock front T does 
not exceed ~ 3 ■ 10*^ K. ISnsman & Burrowsl ( Il992h obtained T 
higher by almost two orders of magnitude, while the color temper- 
ature of the radiation agrees satisfactorily with our value in FigJlSI 
What is the cause of t his contradiction? The poin t is that the de- 
scription of opacity in lEnsman & Burrowsl ( Il992h was too crude. 
These authors simply subtracted the Thomson opacity from the 
Rosseland mean xr and assumed this quantity to be equal to the 
Planck mean for true absorption Qp. Since the Thomson scattering 
at high temperatures dominates in xr and is higher than absorption 
by a huge factor, the subtraction of the two close numbers leads 



to very large errors. In addition, the identificat ion of th e Rosseland 
and P lanck means is an improper procedure dZeldovi ch & RaizeJ 
1 19661) . In our calculations, the true absorption is calculated sepa- 
rately from the total extinction, i.e., always fairly accurately, and 
no problems with the averaging of o pacity over the entir e spec- 
trum arise in multigroup calculations. lEnsman & Burrowsl (Il992l) 
explain the high matter temperature T by the formation of a vis- 
cous jump in transparent layers with strong gas heating (as was 
discussed above, there is no viscous jump in deep layers in a strong, 
radiation-dominated shock). 

To check the influence of opacity on T in a shock, we carried 
out a numerical experiment by replacing the complete opacity cal- 
culation algorithm with the following approximate procedure. The 
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Figure 13. Velocity profiles as a function of Lagrangean mass (top) 
and optical depth r (bottom) in the model 14E1 for SN 1987A from 
iBlinnikov et all fcOOOl) . The time in seconds is given near the curves. 



Figure 14. Matter density for the variant 14E1.3 at the epoch of shock 
breakout versus Lagrangean mass Mr (top) and Eulerian radius r (bottom). 
The proper time is given near the curves. 



code that we developed for the equation of state solves the Saha 
system of equations for an arbitrary number of elements at all ion- 
ization stages. In the calculations that we will now descri be, the 
ionizati on is described in the "mean-ion" approximation by lRaizej 
( Il959ah (see also the book bv lZeldovich & Raizerl ( Il966h ). 

We obtain the continuum opacity Xc needed for our calcula- 
tions using the following procedure (only in this numerical experi- 
ment!). Initially, we take an expression f or the monochr omatic ab- 
sorption coefficient similar to that used bv lVitenselligSlh and based 
on the Kramers approximation for our mixture of elements. This 
expression, plus the Thomson opacity, is not used directly as Xc but 
is used only to obtain the Rosseland mean xr through approximat e 
integration proposed by ( lRaizeij[l959bt Izeldovich & Raizeijll96d) . 



W e obtain the true absorption coefficient ct-p just as was done 
by lEnsman & B urrows ( 199_2), thereby qualitatively reproducing 
the algorithm from their paper, and use Qp instead of the coixect 
monochromatic absorption coefficient in our multigroup calcula- 
tions. We will call the opacity obtained in this way Zeldovich- 
Raizer one and will denote it by ZR. 

The matter velocity and temperature for the variant 14E1.3 
with ZR opacity at the epoch of shock breakout are shown in 
FigJ16l We immediately see that the matter in this variant speeds up 
to higher velocities and heats up in the viscous jump to enormous 
temperatures, T ~ 10® K. Both these effects are explained by un- 
derestimated true absorption: the quasi-adiabatic acceleration of the 
matter in a self-similar shock dCandelman & Frank-Kamenetskvl 
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Figure 15. Bolometric light curve (in the inset) and color (solid line) and 
effective (dashed line) temperatures for two versions of the model 14E1 (see 
the text). 
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Figure 16. Matter velocity (10"^ km/s, top) and temperature (bottom) for the 
variant 14E1.3 with Zeldovich-Raizer (ZR) opacity at the epoch of shock 
breakout versus Lagrangean mass A/^. The proper time is given near the 
curves. 



ll95dlSakura"illl960l) continues longer (cf. the velocity profile in 
Fig J lib, while the matter behind the density jump cannot efficiently 
give up heat and heats up virtually to values corresponding to an 
adiabatic shock. 

Figure [17] compares the calculated bolometric light curves 
for model 14ElshbaX radaNhm23 of the SN I987A presupemova 
dBlinnikov et alj200(]|) . The contribution from the RADA algorithm 
to the computation manifests itself mainly in a careful allowance 
for the radiation time delay, which leads to a decrease in radiation 
flux at maximum and to a broadening of the light curve peak. 

It should be noted that light curve calculated by RADA algo- 
rithm does not have second peak presented earlier dTolstovlbOlOl) 
due to more accurate allowance for boundary conditions. 



4.3 The dependence of the parameters of the flash from the 
explosion energy and the radius of presupernova 

In the models for SN 1987 A, the explosion energy was varied at 
constant presupemova radius. In addition, we constructed two more 
models with greatly differing radii, 10 and 300 Rq. The presuper- 
nova model 14E1 for SN 1987 A with initial radius Ro = AS Rq, 
ejection mass M — 14.7 Mq, and iSkin = 1-0 ■ 10^^ erg was taken 
as the initial one. The explosion energy in these models was the 
same. We constructed models with radii slightly different from the 
standard value: decreased i?o = 40 i?Q(e.g., 14E1.34R) and in- 
creased i?o ~ 58 i?Q, model 14E1.4R6. In addition, we also con- 
structed models with a great increase in the radius to 300 J?0and 
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Table 3. Predictions for the first brightness peak 



variant 


t,d 


ibol> 


erg/s 


Tc,K 




K 




2/3. cm 


r^'^ Ldt, erg 
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14E0.7 
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1012 
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14E1 
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.06726 
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1.339 106 
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14E1.34R 


.05620 


9.175 
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10^ 


2.74 
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10"^ 
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5.30 


10^ 


3.95 
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Note. The T^g maximum almost coincides in time with the Lboi peak, while the Tc maximum occurs ~ 100 s 
earlier 




7400 



Figure 17. Comparison of the bolometric light curves at the epoch of shock 
breakout for an SN1987A-type presupemova model. The solid and dotted 
lines represent, respectively, the STELLA and RADA calculations in the co- 
moving frame of reference. Dashed line represents RADA calculations in 
observer's frame of reference taking into account radiation time delay. 



with its decrease to 10 Rq. The density distribution in the new 
models was obtained homologically from the model 14E1. The out- 
burst parameters are presented in Table. |3]and|4] 

The kinetic energies of the ejecta with radii of 10 and 300 Rq 
slightly differ, because to achieve complete coincidence between 
their values, the calculation should be performed several times by 
iterating explosions with trial energy release at the center. We see 
from FigJlSlthat the outburst time naturally increases with radius 
(for a discussion of this dependence, Blinnikov et al, 2000),while 
the peak luminosity depends weakly on radius: when the radius in- 
creased by a factor of 30, from 10 to 300 Rq the luminosity peak 
decreased approximately by a factor of 2. The effective tempera- 
ture maxima for three versions of the model 14E1 with different 
radii shown in FigJlSldepend more strongly on Rt^: in accordance 
with the definition of L = AnoT^ff^RQ, we have T^s cc Rq^^^ 
at constant L. Recall that the color temperature is approximately a 
factor of 3 higher than the effective one in these outbursts, as can 
be seen from table [3] 

An important parameter of a supernova is the maximum ejecta 
velocity. For various versions of the series 14E this parameter is 




0.2 0.3 
t, days 
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Figure 18. Absolute bolometric magnitudes (top) and effective temperature 
(bottom) for three versions of the model 14E1 with different radii at the 
epoch of shock breakout. The radii are given in solar units. 
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Table 4. Maximum matter velocity for various versions of the series 1 4E at 
an ejecta mass of 14.7 Mq 



Bkin, 10" erg 


Ro, Rq 


"max, lO^km/s 


1.6 


48 


42 


1.3 


48 


37 


1.0 


48 


33 


0.7 


48 


28 


1.2 


10 
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48 


37 


1.1 


300 


20 



Table 5. Maximum matter velocity in various approximations for opacity 
in the model 14E1.3 at an ejecta mass of 14.7 Mq 



Opacity 


Ekin, 10" erg 


Ro, Rq 


limax, lO^km/s 


Standard 


1.3 


48 


37.2 


'high' 


1.3 


48 


37.7 


absorption 


1.3 


48 


37.7 


ZR 


1.3 


48 


60 



the coherent Thomson scattering, while the photon energy could 
change only through the d ivergence of bulk motion, just as in 
iBlandford & Pavnel I198T3) through the term; 



(11) 



where v is the photon frequency, u is the matter velocity, and / is 
the mean photon occupation number (to be more precise, the zero 
angular moment of the occupation number). In all of the remaining 
versions, this term, of course, was also present in the equation for 
the photon energy density. 

Le t us compare in details the equations of lSlandford & Pavnd 
( Il981bl) with equations used in our calculations. 

The monochromatic radiation energy equation and monochro- 
matic radi ation momentum equation with the accuracy 0{v/c) is 
written as jMihalas & Mihalaslll984 eq. 95.18-95.19); 



DEo{vq) 
Dt 
Ov 

+ — [Eo(vo) + Po(i^o)] + -^Foii^o) ■ 
or 



2a . 



d 
duo 



uo(-[Eo{iyo) - Po(i'o)] + ■T^Poi'^o) + 
\r or 

= 4nrio{uo) - cx()(«^o)-Eo(^'o) 



(12) 



given in Table. |4] Table. |5] shows the influence of various assump- 
tions about opacity on this parameter. Standard opacity suggests al- 
lowance for spectral lines in the approximation of a static medium 
(i.e., the values that correspond to the pre-shock region are taken). 
"High" opacity is the case where the spectral lines are everywhere 
broadened by the maximum velocity gradient. Absorption is the 
case where the entire extinction is treated as true absorption. We 
see that there is almost no difference in maximum velocity in all 
these versions. This is because the contribution of lines to the total 
extinction in a hot medium is small. A great difference is obtained 
only for ZR opacity, which underestimates the true absorption on 
the photoeffect. According to the Kirchhoff law, weaker absorption 
also means weaker emission, a case closer to the adiabatic selfsim- 
ilar one, whence the high ejecta velocity. 



4.4 The influence of a frequency grid and approximations in 
treating the Compton effect 

It is particularly interesting to trace whether high temperatures can 
be reached when more realistic opacities than those in the ZR case 
are taken into account for a more careful allowance for the photon 
production, absorption, and scattering processes. In particular, the 
cooling of electrons through the inverse Compton effect (see, e.g. 
IZeldovichI ( [l975l) ) should be taken into account; 



Wcs 



nleC 



(10) 



The results presented here for the models 14E took this ef- 
fect into account roughly — in the approximation of complete 
photon thermalization even at the first scattering, i.e., the photons 
were assumed to be absorbed with a cross section OTk-sTe / rrieC^ . 
To check whether this treatment leads to an underestimation of 
the temperature, we carried out the following experiment. Among 
the models for SN 1987 A, we performed a series of calcula- 
tions for the 14E1X2- type versions, where such thermalization 
of photons was switched off. Here, we took into account only 



Dt 
' dv 



dr 



+ 



C2 
d 

dvo 



(V 1 Dv 



+ -Po{vo, 



= Xo(fo)-Po(fo) 

c 



(13) 



Here Eo{uo), Fo{vo), Po(t'o), Po(j^o) - moments of the radiation 
field, the affix "0" denotes the comoving frame, v - the fluid veloc- 
ity, X - the extinction coefficient, rj - the emission coefficient, the 
time derivatives (D/Dt) are evaluated in a moving fluid element. 

STELLA solves the equations, which are equivalent to equa- 
tions 1 12ll3t, with the neglect of several terms dMihalas & MihalasI 
ll984lcastorll972h ; 



dt 

H — (SlCi, — Ju 
r 



c d 2 
r or 



dr 



{r^u){J^+K,^) 



— — Ju) 2 



d 



(r u)K.v 



(14) 



dt 

-2 



dKLu 
dr 



du 
dr 



(15) 



where u - the fluid velocity, Xa = Xa(p, T, v) - true absorption co- 
efficient, Xs - scattering coefficient, Jv^Hv^tCv - moments of the 
photon occupation number and, for example, Ei, = j (? Jy, 

^"diit ~ numerical stability term. 

In the diffusion limit when K, — >■ j7/3 at photon mean free 
path Ap — the equations I ll4ll5t lead to: 



djy 
dt 



V 1 
3 r2 



d , 2 ^dJv 
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d_ 

dr 



+ C(7?. -XaJ.), (16) 



3(Xa + Xs) dr 

which is equivalent to equation derived by iMihalas & Mihalai 
( 1 1984 eq. 97.68) for the diffusion limit: 
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dEo{uo) D (I 



= V 



+ 



(17) 



di/Q Dt\p)\ V^xoivo) 
+Ko[vo)[A-nB[va,T) - cEo[vo)] 

The left and right part of the equation ( 116) correspond to the left and 
right part of the equation ( 117) and it is supposed that Xv = 
and riv = fc^_B^ + a^Jv Here J,j is mean intensity or zeroth mo- 
ment of radiation field over angles and Ei, = (47r/c)Jy, Bi, - 
isotropic specific intensity in thermal equilibrium, 77 - emission co- 
efficient, X - extinction coefficient, ki, - true absorption coefficient, 
(T - scattering coefficient, p - density of the matter, and recalling the 
equation of continuity D / Dt = d / dt -i- vd / dr 

{Dlnp/Dt) = -r'^[d{r'^v)/dr] = -{dv)/idr) - (2i;/r).(18) 

The equation used bv lBlandford & Pavnd ( fl981bl , eq.l8) can 
be rewritten as follows (we omit the photon energy redistribution by 
Compton scattering due to STELLA does not take it into account): 



dn 



u ■ Vn - 



1 ,„ . dn 



Vn + j, (19) 



dt 3 dv \Zneo{u) 

where n - mean occupation number, i.e. J in the equation ( 116) 

The first two terms of the equation ( 119) correspond to the first 
term of equation ( 116) (the time derivative in STELLA is taken at 
fixed Lagrangean radius). The third term of (119) corresponds to the 
second teim of equation ( 116) and describes the heating of the radia- 
tion through compression. The first term of the equation ( 119) in the 
right part corresponds to that one in the equation ( 116) and describes 
diffusion of photons in space. The photon source term j in equation 
( 119) is added without without any corresp onding absorption term 
and it is inconsistent with thermodynamics dPsaltis & Lambll997h . 
The equation J16) contains both emission and absorption terms: 
c(?7i' - X^Jy)- 

We note that our equations in clude all terms of final equa- 
tion of iBlandford & Pavnd ( Il981bl) , but our approach is more ac- 
curate because_^U_varmWes_ar^ in Lagrangean frame, 
while Isiandford & Pavnd ( Il981b|) did not dist inguish clearly the 
fluid and the inertial frame ( Fukue et al.lll985l) . Moreover, we do 
not add any inconsistent terms and do not confine ourselves to the 
diffusion limit. 

The characteristic time scales of the changes in spectrum due 
to Compton scattering f^^, f^^, and the divergence of bulk motion 
are, respectively. 
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(20) 
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(22) 



where the scale length ~ c/{neC7Tu) dBlandford & Pavn3ll981bl) . 
In the case of ultrarelativistic motions (see the transfer equation ([4J 
below in the comoving frame of reference), the term with udf /dv 
takes a more complex form (the intensity is I — 2hv^ f / c?). 

Another difference of the 14ElX2-type versions is a wider fre- 
quency grid: the number of frequency bins was doubled and the 
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Figure 19. Spectral flux distributions uL^, for four instants of time (shown 
in seconds) for the variant 14E1X2 at shock breakout (solid line) and the 
best fit by a blackbody spectrum(dashed line). The time lag was disregarded. 



minimum wavelength A was set equal to A was set equal to 10^^ A 
(instead of 1 A in the standard series 14E). 

The resulting outburst light curves for two versions of our cal- 
culation are presented in Fig. [15] We see that the difference is very 
small. Fig.[T9]shows typical spectra at the epoch of shock breakout 
for 14E1X2. Increasing the grid shows that the shock breakout in 
X rays becomes visible much earlier than in visible light, because 
the absorption is weaker there. However, the fluxes in the hard part 
of the spectrum remain low. The figure also shows the best fit by a 
blackbody spectrum with a temperature that we call the color one. 



4.5 The hardest semirelativistic variant, the SN Ib/c model 

As we see from Table |4] the maximum ejecta velocity depends 
strongly on presupernova radius: it is higher in more compact stars. 
Therefore, it is interesting to consider the shock breakout for type 
Ib/c supemovae (their progenitors are Wolf-Rayet stars with radii 
of the order of the solar one or even smaller). The type-Ib/c pre- 
supernova model that we use was constructed with the KEPLER 
code by an evolutionary computati on from a main-sequence star in 
IWooslev. Langer & Weaveil ( ll995h . the model 7^4. At the end of its 
evolution, the presupernova star has a core composed of helium and 
heavy elements with amass of 3.199AfQ and a radius of 1.41 • 10^^ 
cm. The radius in this case was fixed "manually", because in the 
outer stellar layers KEPLER models the stellar wind and the model 
is not in hydrostatic equilibrium. 

This model was chosen because the velocity of matter at shock 
breakout reaches the values about a half of speed of light. At such 
velocities, the relativistic corrections to the radiative transfer are 
significant and applying the RADA code here is of great interest. 
Although STELLA does not work with hydrodynamic flows with 
large Lorentz factors, it takes into account the relativistic effects 
with an accuracy of iv/c) sufficient for the models we consider. 
The results of calculations of light curves and spectra of this model. 
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Figure 20. Comparison of the bolometric light curves at the epoch of 
shock breakout for a type Ib/c presupemova models: explosion energy 
E = 9 ■ 10^^ erg, maximum matter velocity Vmax ~ 0.3c (top) and 
E = 1.8 ■ 10^^ erg, Vmax ~ 0.5c (bottom). The solid and dotted lines 
represent, respectively, the STELLA and RADA calculations in comoving 
frame of reference. Dashed line represents RADA calculations in observer's 
frame of reference taking into account radiation time delay in the observer's 
frame of reference. 



which we denote slb7a2, in the variant of the algorithm RADA with 
a soft spectrum (A > lA) are the following. 

Figure l20lcompare the calculated light curves for the model. 

We see from the computational data that allowance for the de- 
lay effect and a strict allowance for the relativistic radiative transfer 
affect the light curve shape. It leads to decrease in radiation flux at 
maximum and to a broadening of the light curve peak. The effect of 
geometrically eclipse the radiation of the outburst from the edge of 
the star (see Fig. lB2t i in Appendix B) is revealed in the calculations 



Figure 21. Matter temperature for the variant slb7a2X at shock breakout 
versus Lagrangean mass A/,, measured from the surface. The time in sec- 
onds is given near the curves. The temperature peak is at an optical depth 
T~ 200, 50, 4, 1, 0attimes 67.0, 67.5, 67.9, 68.1,68.9 s. 



but it is quite small in this model due to small velocities after shock 
breakout (0.2-0.5 c) and can not be seen on the graph in details. 
The light curve pres ented here is c alculated more accurately than 
we discussed earlier ( lTolstovll2010l) . 

Below we consider slb7a2 variants of STELLA with a hard 
spectrum (A < lA) as well. Options differ from that described 
above, only by the numerical treatment of radiation transfer. 

Program STELLA, which is based on the nonrelativistic equa- 
tions of hydrodynamics, can be modified to allow for relativistic ve- 
locities. Instead of the velocity u, the numerical scheme may work 
with the quantity U = 711. The equati ons for U appear a lmost as 
in the Newtonian limit for u (see, e.g., IUrzhumovl ( l2002h ) but the 
velocity u = U /[U"^ /c? + lY^'^ now can not exceed the speed of 
light. In the current work we do not use this modification and limit 
ourselves to the case when u < c/2. 

In all the versions of the model slb7a2, described above, we 
changed neither the explosion energy nor the radius. The variant 
slb7a2X was computed by STELLA on the same spatial grid as 
slb7a2, but the number of frequency bins was doubled and the 
minimum wavelength A was set equal to 10^"^ A (instead of 1 A 
in slb7a2). In this case, no photon thermalization under the Comp- 
ton effect was imposed. 

In Fig. [2TI the matter temperature for the variant slb7a2X is 
plotted against the Lagrangean mass Mr measured from the sur- 
face. We see that the maximum temperatures are enormous — up 
to 10^'^ K (i.e., of the order of MeV). The spectrum for an observer 
in the comoving frame at the surface is shown in Fig.|22] Since the 
STELLA algorithm includes the evolution of photons in a converg- 
in g flow in the shock in the sam e approximation dl It . as considered 
by blandford & Pavnd [l981ij) . one could think that our computa- 
tion confirms their analytical result. 

In fact, however, the high density and hardness of the radiation 
in such a shock obtained in the slb7a2X computation require a 
more careful allowance for the photon production, absorption, and 
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Figure 22. Spectral flux distributions vLv for four instants of time (shown 
in seconds) for the variant slb7a2X at shock breakout (solid line) and the 
best fit by a blackbody spectrum(dashed line). The time lag was disregarded. 



Figure 23. Matter temperature for the variant slb7a2Xm6 at shock break- 
out versus Lagrangean mass Mr measured from the surface. The time 
in seconds is given near the curves. The temperature peak is reached at 
T ~ 200 50, 4, 0.5 at times 67.0, 67.5, 67.9, 68.1 s, respectively, and 
virtually disappears subsequently. 



scattering processes. In particular, the cooling of electrons through 
the Compton effect, as was discussed above for the variant 14E1X2 
as an example, should be taken into account. 

Nevertheless, it turn s out that allowance even for the weaker 
double Compto n effect jMandl & Skvrmd 1 19521 : IWeaved 1 19761: 
lLightmaij[l98ll) leads to a sharp change of the results. We will 
postpone its rigorous allowance until a future paper, while here this 
effect was simulated by a very small admixture of true "gray" (i.e., 
frequency-independent) absorption, 10^^ of the Thomson scatter- 
ing in an SN lb progenitor. We called this variant slb7a2Xm6. 

The plots in Figurel23landl24lshow that the matter temperature 
and spectrum hardness for the variant slb7a2Xm6 decrease sharply 
compared to slb7a2X. 

The previous results were presented for an observer at the 
"edge" of the supernova ejecta in the comoving frame. The STELLA 
results can be carefully transformed to the rest frame using the 
RAD A algorithm. Here, Fig.|25]compares the monochromatic light 
curves at shock breakout in the reference frame of an observer 
at rest for a type Ib/c presupernova. The gray line represents the 
STELLA computation without the time lag. The black line repre- 
sents the computation in the observer's frame of reference with a 
full allowance for the time lag (in the blackbody approximation for 
brightness) computed by the RADA algorithm. 

Note that the flux in the computation including the time lag 
was calculated under the assumption of blackbody intensity. This 
assumption allows the aberration and Doppler effects to be prop- 
erly taken into account, but it leads to a difference in fluxes before 
the shock breakout. On the other hand, the calculation of fluxes by 
disregarding the time lag is based only on the work with the radia- 
tion times. Therefore, the comparison in the figure may be consid- 
ered only as a qualitative one, demonstrating the main differences 
and peculiarities of the light-curve shapes. 



5 CONCLUSIONS 

In this paper, we found that the high-temperature peak behind the 
shock front and the possibility of the formation of a hard power-law 
"tail" in the radiation spectrum are suppressed through a very weak 
process of photon absoiption and production with a cross section 
at a level of one millionth of the Thomson scattering cross section. 
At this level, the absorption can be provided by double Compton 
effect. 

In theoretical models of type Ib/c supemovae shock breakout 
the motion of matter at velocities near the speed of light should be 
taken into account. Development of the numerical algorithm RADA 
for solving radiative transfer with the hydrodynamic part of the al- 
gorithm STELLA provides a reasonable allowance for the relativis- 
tic effects. 

Our calculations of supernova shock breakouts can be used 
for evaluating and interpreting the detection of supernova explo- 
sions in planned space experiments. In the near future the number 
of observed outbursts will grow due to the launches of new space- 
crafts and the theoretical models are called for because they can be 
used to predict the number and nature of observed outbursts. 

One of these exp eriments could be LOBSTER 
jCalzavara & Matzned l2004t) . which is seems to be discarded 
now. Anyway, the characteristics of a detector of this type are 
more favorable for the detection of outbursts from presupernovae 
that are red and blue giants. For Wolf-Rayet stars, the detection 
rate is estimated to be several outbursts per year from a maximum 
distance of 50 Mpc. This distance is a natural limit of the detection 
of outbursts of this kind, because the daily detection limit of 
the spacecraft detector would be about 10^^^ erg cm~^s~^ 
jPriedhorskv, Peele&Nug"en3ll996l) . For events in our Galaxy, the 
outburst spectrum will be quite resolvable and our computation of 
model 7 A, already refined by relativistic effects, can be used for 
comparison with the observational data. 
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Figure 24. Spectral flux distributions uLv for four instants of time (shown 
in seconds) for the variant slb7a2Xm6 at shock breakout (solid line) and 
the best fit by a blackbody spectrum (dashed line). The time lag was disre- 
garded. 



One of the most interesting problems of the stellar evolution 
theory is an unknown number of outbursts from collapsing SN 
1987A-type supemovae. They are difficult to observe due to the 
short outburst duration (~ 100 s), but the outburst brightness makes 
it possible to detect them with X-ray detectors. 

Type Ib/c supernova explosions are also difficult to observe 
due to their short duration (~ 10 s), but, in this case, the radiation 
has the hardest spectrum and the X-ray radiation from the shock 
propagating within the presupernova wind can be detected. The de- 
tection of type Ib/c explosions is very important for the theory of 
hypemovae, since it will allow us to describe better the connection 
of supernovae with gamma-ray bursts and, in the case of nearby 
supemovae, there can be a correlation between the events and the 
detection of gravitational waves and neutrinos. 

According to the estimations of Calzavara & Matzneil (|2004|) . 
LOBSTER would have the following possibility of outbreak reg- 
istration for 3 years: from 50 to 600 supemovae type II and some 
type Ib/c within 250 Mpc. These estimations are based only on 
several existing calculations taking into account the multigroup ra- 
diation transport descr ibed in the litera ture - the model SN 1 993 J 
l lBlinnikovetai]|l998h and SN 1987 A jBlinnikov et"ai]l200(]|) , and 
the rest of the spectmm of data obtained by scaling these calcula- 
tions. To improve them, we can use the theoretical modeling meth- 
ods we developed. 

A further development of numerical one-dimensional super- 
nova explosion models requires both a more accurate allowance 
for the hydrodynamic effects and a more complete description 
of Compton scattering of radiation. The observed asymmetry ef- 
fects require developing multidimensional algorithms and radiative 
transfer calculations, but the success achieved in modeling super- 
nova explosions serves as an impetus for future studies. 
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Figure 25. Comparison of the monochromatic light curves at shock break- 
out for a type Ib/c presupernova. The gray line represents the STELLA com- 
putation in the comoving frame of reference with the only correction for the 
observer's time t t — R{t) /c. The black line represents the computation 
in the observer's frame of reference with a full allowance for the time lag 
(in the blackbody approximation for brightness). 
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APPENDIX A: TIME-DELAY-SPREAD LIGHT CURVES 

The light radiated by different places of stellar photosphere at the 
same time t reaches observer at different times spread over time 
interval of width R/c, R and c being the photosphere radius and 
speed of light, respectively. The observer actually sees a convolved 
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t R/c 

Figure Al. Time-delay-spread (5-function. 



signal radiated by photosphere within the time interval At — R/c. 
Mathematically, the observed luminosity L{t) can be expressed 
through the intrinsic, not spread, luminosity Lq {t) by the following 
integral: 



Lit) 



2 / Lo{t')xdx, 
'o 



t-{R/c)-{l-x). 



(Al) 



This expression is based on 3 basic assumptions. First, it implies the 
distance D between observer and the star is very large: R/D <C 1. 
Second, it assumes that the radiation intensity I is isotropic all over 
the photosphere — i.e., I does not depend on the angle 9 between 
the direction of radiation propagation and the perpendicular to the 
photosphere surface. Third, the Eq. (lAU neglects the change in R 
during the time interval At = R/c. This is a good approximation 
as long as the velocity of photosphere, V = dR/dt (do not con- 
fuse with that of matter crossing the photosphere!) remains small 
in comparison with the speed of light, say \V/c\ < 0.1. This is true 
for the shock wave breakout in common supemovae of Types 11, lb, 
and Ic. 

Integrating both sides of Eq. dAlt by t, one can make certain 
that L{t)dt 

— Lo(i)dt- Thus as it should be, the transformation in question 
conserves the total radiated energy — it only turns out to be redis- 
tributed in time for a remote observer In a trivial limit i? — >■ (no 
time delay at all), Eq. dAlt gives, naturally, L{t) = Lo{t). The fol- 
lowing two simple but instructive examples illustrate how Eq. l lAlt 
actually works. 



Useful Examples 

Consider the case when the intrinsic luminosity time-scale r is very 
short as compared to i?/c, or a = R/(ct) — oo. In such a limit 
Loit) can be represented as the (5-function: io(t) = E 5{t), where 
E is the total radiated energy. Inserting this La in Eq. l lAlb we ob- 
tain 



L{t) = 



0, 
c 
Ti 

0, 



t 0, 

t> ^. 

c 



(A2) 



The resulting convolved luminosity L{t) has a triangle form as 
shown in Fig JAll 

Another example providing the integral in Eq. dAU in a closed 
form is the exponential function: 




Figure A2. Time-delay-spread exponential function. 

Shown by a dashed curee is the initial Lo(t) (Eq. lA3) that corresponds to 

a = 0. 
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The result is: 
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Figure lA2l shows L(t) for a number of values of parameter a = 

R/{ct). 

The next Fig JA3l shows L(t) in terms of Ec/R as a function of 
time in terms of R/c. The dashed curve represents the (5-function 
response (Fig. lAlt in a limit a — >■ oo. We see that the curve for 
Q = 8 is already not so far from this limit. 

It is easy to find that the time-spread light curve (Eg. lASt 
attains a maximum Lm at t — tm within the time interval 
Q<t< R/c: 

Lm = ^ - [a - ln(l + a)] , = r ln(l + a) . (A5) 

To find the light curve width Atq at a given level L{t) / Lm = q < 
1 one has to specify a and solve a couple of transcendent equations. 
Table l A 1 1 presents Atq for q — 10^"'^ (one stellar magnitude be- 
low maximum) along with Lm and tm for the values of a shown in 
Figs.lA2landlX3l 

The dependence of Atq on a can be approximated with a good 
accuracy (better than ~ 2%) by the following simple formula: 



Ata 



0.921 +0.6 a 1 + 



0.45 

1 + 



, < a < oo . (A6) 



Finally, it is worthy to notice that for time t > R/c the ratio 
A = Lit) / La{t) (also shown in Table lAlt is always greater than 1 
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Table Al. The properties of the time-spread exponential light curve. 



a = R/{ct) 


0.0 


0.2 


0.4 


0.6 


1.0 


2.0 


4.0 


8.0* 


L„,/{E/t) 


1.0 


0.884 


0.794 


0.722 


0.614 


0.451 


0.299 


0.181 


tm/r 


0.0 


0.182 


0.336 


0.470 


0.693 


1.10 


1.61 


2.20 


Atg/r 


0.921 


1.07 


1.22 


1.36 


1.63 


2.29 


3.67 


6.33 


A 


1.0 


1.07 


1.15 


1.23 


1.44 


2.19 


6.20 


92.9 



In the limit a 



oo (see FigjM): im ^ 



■ 2, 



R 



■■ 0, Atg 



: 0.6ra = 0.6 



R 



LlJ 1.5 



1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



a = R/(cT) 
Lo(t) = E 6(t) 




Figure A3. Same as in Fig JA2l but for the luminosity and time in terms of 
Ec/ R and R/c, respectively. 



and does not depend on time. 
L{t) _ 2 



exp(a) 



= 1 + S + ^ + 
^ 3 ^ 12 ^ 



R 

c 



<t < oo. 



(A7) 



Thus, the time-spread luminosity "tail" turns out to be brighter than 
that of the intrinsic luminosity Lo{t) by a factor of A (a). For a > 2 
this effect increases the apparent brightness by about one stellar 
magnitude or more. 

In the case of a power-law light curve, I/o(i) ~ ^1 + , 
that also allows the integral to be taken in a closed form, such an 
effect is absent: L{t) / Loit) — ^ 1 for f — >■ cxj. The exponential 
function decreases fast enough to make the luminosity tail remain 
ever overluminous. 



APPENDIX B: TIME DELAY RELATIVISTIC EFFECTS 

There exist a number of constraints on the lower limit of integration 
/imin in Eq. l|7j and, accordingly, there are several reasons why the 
light from the source does not reach the observer. 

(i) The first constraint is a geometric one determined by the an- 
gular sizes of the radiation source visible to the observer. For ex- 
ample, the back of the surface of the spherical envelope at rest is 
disregarded in the integration, which corresponds to a zero lower 
limit of integration for any time (Fig. IB 1 fa ). 

(ii) The second constraint is a dynamical one related to the 







Figure Bl. Dependence of the radiation flux at a remote point of observa- 
tion (on the right) on the dynamical characteristics of the observed envelope. 
The envelope region from which the radiation reaches the observer is high- 
lighted by the heavy line for the following cases: the envelope is at rest (a), 
expands with a constant velocity v = jSc (b), and expands with accelera- 
tion (c). Panel (d) shows the equitemporal surface the radiation from some 
regions of which does not reach the observer due to the peculiar envelope 
dynamics. 



finiteness of the speed of light. If the radius R{t), where radiation 
begins to freely propagate outwards, does not depend on or shrinks 
with time then /imin = 0. However when R{t) increases with time 
the lower limit of integration becomes ^min = /3o = - dR{t)/dt 

(Fig.mij). 

(iii) An even more severe constraint can be in the case of accel- 
erated envelope expansion. In this case, it can reach some spatial 
point faster than the light emitted from the edge (relative to the 
observer) regions of the envelope (Fig. IBU . The light from the en- 
velope point {R{t),iJ,) does not reaches the observer in case the 
following condition is performed: 



3dt : G{t, n, dt) = 

= R^{t) + [cdtf + 2R{t)p,cdt - 



R^{t + dt)<0. (BI) 



After all these points are determined the flux for the distant ob- 
server can be found by integration over the equitemporal surface 
from which the photons reach the distant observer simultaneously 
and by excluding all found points from the integration (Fig. IB IB ). 
This constraint is not active if second derivative of the function 
G{t, fj., dt) is always positive. 

The last constraint is most stringent and includes the previous 
constraints as special cases. Let us demonstrate this using the fol- 
lowing example. For the dynamics of an expanding envelope with 
a linear increase in velocity, the propagation law can be written as 
R/c — Ro + Pot + at^, where Rq, /3o, and a are constants. Gen- 



© 2012 RAS, MNRAS 0Q0.mi20l 



20 A.G. Tolstov, S.I. Blinnikov, D.K. Nadyozhin 



erally, the analytical solution is cumbersome and cuts out certain 
intervals of angles from the flux integration. However, it is reduced 
to /Xmin = 1^0 + lot j c for i?o = and we can see that, in this case, 
the flux will be lower than that in the case of envelope motion with 
a constant velocity. 

To demonstrate the peculiarities of the flux behavior at the 
point of observation, let us consider the effect of a flux decrease 
in the problem of an outburst of radiation on the envelope surface. 
This problem for an instantaneous outburst and an env elope at rest 
was solved bv llmshennik, Nadezhin & UtrobinI ( Il98ll) . It is simi- 
lar to the pr oblem of a light echo fro m a nova outburst considered 
long a go bv lCouder3 ( Il939h . see also lBarkov & Bisnovatvi-KoganI 
( l2005h for an infrared outburst generated by a gamma-ray burst in 
a dust cloud. In contrast, we investigate the influence of an outburst 
with duration dt followed by envelope expansion on the radiation 
flux at the point of observation; the source cannot be considered 
stationary and instantaneous and the envelope motion can be rela- 
tivistic. 

Consider a spherical envelope with radius i?o = c on the sur- 
face of which an outburst with intensity Iv = 100J^,o and duration 
dl = 0.1 s occurred at t = followed by envelope expansion 
according to the law i? = i?o + 0.5c(t — d£). The calculated ra- 
diation flux from the envelope at a remote point of observation is 
shown in Fig. lB2l (scenario 1). Also shown here is a comparison 
with the flux calculated for the case where only the outburst with- 
out subsequent motion occurred on the envelope surface (scenario 
2) and the flux calculated when there was no outburst but the ex- 
pansion followed (scenario 3). The gentler flux decline in the first 
scenario than that in the second one stems from the fact that the flux 
from the envelope starting to move at f = 0.1 s is higher that from 
the envelope at rest. The steep decline near t — 0.6 s is related to 
the fact that the light from the moving parts of the envelope arrives 
at the point of observation earlier than that from the parts at rest 
and hy t = 0.8 s the moving parts completely obscure the delayed 
flux of the outburst from the envelope edges and the flux behavior 
begins to follow scenario 3. 
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Figure B2. Radiation flux at a remote point of observation from an enve- 
lope of radius Rq = c on the surface of which an outburst with intensity 
Iv = 100/y,o and duration dt = O.IRq/c s occurred at i = followed 
by envelope expansion R = Ro + 0.5c(t — dt) (scenario 1). The dotted 
line denotes the flux as if only the outburst without any subsequent expan- 
sion occuiTed (scenario 2) and the dashed line represents only the motion 
(scenario 3). A zero time of observation corresponds to the outburst onset; 
Fq is the radiation flux at a remote point of observation at the initial time 
t = 0. 
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